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We investigate the nucleosynthesis in a massive star of 70 Mq with solar metallicity in the 
main sequence stage. The helium core mass after hydrogen burning corresponds to 32 Mq. 
f~i ' Nucleosynthesis calculations have been performed during the stellar evolution and the jetlike 

CL(, supernova explosion of a collapsar model. We focus on the production of elements heavier 

than iron group nuclei. Nucleosynthesis calculations have been accomplished consistently 
from hydrostatic to dynamic stages by using large nuclear reaction networks, where the weak 
S-, p-, and r-processes are taken into account. We confirm that s-elements of 60 < A < 90 



~^ , are highly overproduced relative to the solar abundances in the hydrostatic nucleosynthesis. 

During oxygen burning, p-elements of yl > 90 are produced via photodisintegrations of 

seed s-elements. However, the produced p-elements are disintegrated in later stages except 

^S| ' for ^^''Ta. In the explosive nucleosynthesis, elements of 90 < A < 160 are significantly 

^ , overproduced relative to the solar values owing to the r-process, which is very different from 

QO ■ the results of spherical explosion models. Only heavy p-elements (A'^ > 50) are overproduced 

OO ' via the p-process because of the low peak temperatures in the oxygen- and neon-rich layers. 

^^ I Compared with the previous study of r-process nucleosynthesis calculations in the collapsar 

\^ , model of 40 Mq by Fujimoto et al., [S. Fujimoto, M. Hashimoto, K. Kotake and S. Yamada, 

Astrophys. J. 656 (2007), 382; S. Fujimoto, N. Nishimura and M. Hashimoto, Astrophys. J. 

680 (2008), 1350], our jet model cannot contribute to the third peak of the solar r-elements 

and intermediate p-elements, which have been much produced because of the distribution 

of the lowest part of electron fraction in the ejecta. Averaging the overproduction factors 

over the progenitor masses with the use of Salpeter's IMF, we suggest that the 70 Mq star 

could contribute to the solar weak s-elements of 60 < A < 90 and neutron-rich elements of 

, J , 90 < A < 160. We confirm the primary synthesis of light p-elements in the ejected matter 

|r%J . of high peak temperature. The ejected matter has [Sr/Eu] ~ —0.4, which is different from 

H ■ that of a typical r-process-enriched star CS22892-052 ([Sr/Eu] ~ -1). We find that Sr-Y-Zr 

. y .' isotopes are primarily synthesized in the explosive nucleosynthesis in a similar process of the 

primary production of light p-elements, which has been considered as one of the sites of a 

lighter element primary process (LEPP). 
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§1. Introduction 

The origin of elements, particularly those heavier than iron, is still under de- 
bated Since charged particle reactions are difficult to produce those elements inside 
stars because of coulomb barriers, other nucleosynthesis processes, that is, two neu- 
tron capture processes, are required. One is the r (rapid)-process and the other is 
the s (slow)-process.™ In the r (s) -process, neutron captures are faster (slower) 
than beta decays. Since the r-process requires high neutron exposure relative to 
seeds, the r-process favors low electron fraction (1^) and/or relatively high-entropy 
environments ^ 

One of the promising sites of the r-process has been thought to be the neutrino- 
driven wind.'^''^''^''^ However, recent one-dimensional hydrodynamical simulations 
of the neutrino-driven wind with Boltzmann neutrino transport have revealed^ ' that 
the electron fraction of the wind becomes high (YJ, > 0.5) and the entropy becomes 
low for the r-process. Therefore, other astrophysical sites such as neutron star merg- 
gj.^ioji,[TTj,[T2ji Qj, black hole winds^SJ have been proposed. However, the properties of 
ejecta such as densities, temperatures, and electron fractions are highly uncertain ex- 
cept for those in Ref. I12p . which makes even the qualitative analysis of the r-process 
difficult. 

In general, the s-process occurs at the end of core helium burning in massive stars 
and/or in AGB stars. Elements heavier than iron of A < 90 are produced in massive 
stars, which is called the weak component of the s-process (weak s-process) .l^ On 
the other hand, elements of 90 < A < 208 are produced in AGB stars called the 
main component (main s-process) i^ The s-process is very sensitive to cross sections 
of neutron captures and /3-decay rates, especially at the branching points such as 
^^Se and '^^Kr.-'-^' Therefore, the nucleosynthesis in massive stars with the use of 
recent experimental cross sections is worth investigating. From the view point of 
astrophysics, part of the elements synthesized by the weak s-process could be ejected 
through the subsequent supernova explosion, where the produced s-elements should 
be the seeds of the p-process.'^ 

Abundances of metal-poor stars provide a good opportunity for understanding 
the nucleosynthesis because the abundances reflect the outcome only from a small 
number of supernova explosions. Observations of metal-poor stars have stron gly 
suggestecpSi that r-elements of the extremely metal-poor stars that have [Eu/Fe]*^l 
> 1 (hereafter referred to as r-process-rich stars) have a "universal" abundance 
pattern, which reproduces the pattern of the solar system r-process abundances for 
Z > 56. However, the abundances of Z < 56 are not the case (e.g., Ref. [T7|) ). and 
the observed [Sr, Y, Zr/Ba, Eu] ratios have dispersion in low-metallicity stars.l^ In 
particular, on the basis of the models of the chemical evolution of galaxies, it has 
been suggested-^J that the abundances of Sr-Y-Zr {Z = 38, 39, and 40, respectively) 
estimated from the contributions of the s- and r-processes are about 10 to 20% less 
than the solar system abundances. As a consequence, a primary component from 
massive stars is needed to explain 8% of the solar abundance of Sr and 18% of 



*■• We adopt the usual notation [A/B] = log(A'A/AO^) — log(A''A/AfA)0 for elements A and B. 
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those of Y and Zr, which should require a so-called lighter element primary process 
(LEPP)P 

On the other hand, the mechanism of core-collapse supernova explosions is 
still a topic of debate. Pushed by recent observations revealing the aspherical na- 
tures of sup er novae j™''^ multidimensional studies of core-collapse supernova ex- 
plosions have been elaborately performed as described in reviews .'^''^''^''^''^ 
Recent two/three-dimensional neutrino-radiation hydrodynamic simulations have 
shown successful supernova explosions, although in some of the models, the explo- 
sion energies are relatively small (10^^-10^*^ gj-g)|26li,|27li,[28li:|2lli qqj) phase transition 
with a mixed phase of quarks and hadrons has also been reported as another pos- 
sible supernova explosion mechanism even though the explosion is assumed to be 
spherical.'^ Magnetohydrodynamical (MHD) simulations with some approximate 
neutrino transport schemes have showEpli''l22}>l331i.lS4),35) jg^jji^g explosions under some 

specific combinations of initial parameters for a strong magnetic field and differen- 
tially rapid rotation. While neutron stars are expected to be left after supernova 
explosions, it has been suggested that a star of more than 25 Mq may collapse 
to a black hole (BH);^^' an accretion disk is formed around the BH if the star 
has enough angular momentum before the collapse. This system could produce a 
relativistic jet of gamma-ray bursts (GRBs^^ZJ due to MHD effects and/or neutrino 
heating around the rotational axis,'^''^ whose system is called a collapsar model!^'''' 
MHD simulations in the context of the collapsar model have shown the formation of 
jei^mmmmmmm due to winding-up effects of the magnetic field or the 
Blandford-Znajek process. ■ 

Nucleosynthesis calculations of the r-process with the use of a collapsar model 
of 40 Mq have been performed extensively by Fujimoto et al.,'^''^ where it is shown 
that the r-process would operate inside the jets. Explosive nucleosynthesis in GRB 
jets has also been investigated .'^''22' Recent nucleosynthesis calculations in a three- 
dimensional MHD supernova model have suggested that such supernovae could be 
the sources of the r-process elements in the early Galaxy.^^' However, in those cal- 
culations, the produced nuclei are limited to primary synthesized ones inside the jets 
and comparisons with the solar system abundances have been focused on elements 
heavier than iron group nuclei. Nucleosynthesis calculations in spherical supernova 
explosions and detailed hydrostatic ones of the progenitors have proved that elements 
of 20 < j4 < 90 are co-overproduced relative to the solar system abundances.'^ How- 
ever, elements of ^ > 90 are not overproduced except for some p-elements. Recently, 
explosive nucleosynthesis calculations for a 15 Mq presupernova model^ with the 
solar metallicity based on two-dimensional hydrodynamical simulations have been 
performed,'^ in which neutrino-driven explosions are triggered by adjusting the 
core neutrino luminosity parametrically. They have concluded that the overproduc- 
tions relative to the solar abundances are similar to the results of spherical explosion 
models.'^ In our previous paper (Paper I),'^ we performed explosive nucleosynthe- 
sis calculations inside the jetlike explosions for the collapsar of a 70 Mq star with 
the solar metallicity. These calculations include hydrostatic nucleosynthesis using a 
nuclear reaction network, which has 464 nuclei (up to ^^Kr). In the present paper, 
we revisit the nucleosynthesis inside the jetlike explosion of the collapsar model and 
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hydrostatic one taking into account all of the weak s-, r-, and p-processes. This 
makes it possible for us to estimate the consistent abundances of the ejecta. In 
particular, we study whether the collapsar model could be the source of elements 
heavier than iron. 

In §2, we present the hydrostatic nucleosynthesis during the evolution from the 
helium burning stage to the onset of the core collapse. In §3, we briefly summarize 
the MHD explosion model for the explosive nucleosynthesis and show the results. 
Section 4 is devoted to a summary of the overall results. We give some discussions 
with respect to uncertainties and remarks about the connection between metal-poor 
stars and the possibility for LEPP. 

§2. Hydrostatic nucleosynthesis 

In Paper I, we investigated the nucleosynthesis in a massive star of 32 Mq 
helium core corresponding to a main sequence star of 70 MqW^ We have used the 
evolutional tracks with a nuclear reaction network, which includes 464 nuclei (up to 
94j^j,^156ji jj-^ massive stars, the (weak) s-process should occur at the end of helium 
and carbon burning stages. However, in the previous paper, we have included only 
nuclei oi A < 94; we could not discuss the weak s-process. Therefore, in the present 
paper, we perform a more detailed nucleosynthesis calculation with a larger nuclear 
reaction network. 

2.1. Stellar model, initial compositions, and physical inputs 

A star of Mms ~ 70 Mq with the solar metallicity could correspond to the 
upper limit of accreting BH models (collapsars), because more massive stars suffer 
from the strong mass loss. ■ As a result, the size of the helium core will be affected 
considerably by the mass loss. However, the rate of mass loss is still very uncertainp^ 
we calculate the evolution of a massive helium core, Mq = 32 Mq, without the mass 
loss as an extreme case, which is worth studying to see the final fate for the series 
of helium core evolution. 

We calculate the nucleosynthesis along each evolutional track of the Lagrange 
mass from the stage of gravitational contraction of the core to the initiation of iron 
core collapse. The calculation has been carried out by using changes in the den- 
sity (/)), temperature (T) and convective regions. This is the so-called postprocess 
nucleosynthesis calculation. In convective regions, elements are mixed and composi- 
tions become almost uniform. Therefore, the region is calculated as one zone with 
averaged mass fractions and nuclear reaction rates as in Ref. I59p . 

Toward s-process calculation, we construct a new reaction network including 
1714 nuclei up to ^^^U, in which the reaction rates are based on a new REACLIB 
compilation, namely, JINA REACLIB database.l^ The included elements are given 
in Table [H The experimental (n, 7) reaction rates in JINA REACLIB are based on 
KADoNiS projects."^ In stellar environments, /3-decay rates could be different from 
the values in laboratories. Takahashi and Yokoi (hereafter referred to as TY87)^ 
calculated theoretically the /3-decay and electron capture rates for elements heavier 
than ^^Ni and tabulated the rates taking into account thermally enhanced ionized and 
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excited states in stellar interiors. The table ranges over 5x10^ < T < 5x10^ K and 
10^^ < Ue < 3x10^^ cm~^, where Ue is the number density of electrons. We adopt 
the temperature and density dependence of the rates if available. Unfortunately, 
the range of the table is limited only for T and p of the helium burning stage. 
If the temperature and density are outside of the table, we use the values of the 
edges of the table. Reaction rates concerning ^^'^Ta are specially treated as noted in 
Appendix A, because ^^"^Ta has a long-lived isomeric state. As suggested by Prantzos 
et al.,!^ the timescales of neutron- induced reactions (w 10~^ s) are much shorter 
than those of the variation of abundances of the other elements (~ 10^-10^^ s) as 
well as convective timescales, and thereby, neutron abundance is determined locally 
under thermal equilibrium conditions. This indicates that neutrons are not mixed 
uniformly in convective regions. Therefore, to calculate abundances in a convective 
region as one zone including the effects of different neutron abundances over the 
region, we adopt doubly averaged reaction rates for "neutron-induced" reactions of 
(n,7), {n,p), and {n,a), according to the same method described in Ref . [59jl . 

Let X{i) denote the mass fraction of the element i. The initial mass fractions 
are assumed to be XC^Re) = 0.981 and X(i'^N) = 0.0137, where all the original CNO 
isotopes are assumed to be converted to ^'^N during the core hydrogen burning. Mass 
fractions of the heavier elements are taken to be proportional to the solar system 
abundanceP* (e.g., X{^^Fe) = 1.17 x lO'^). 

The evolutionary changes in composition with respect to the stellar structure 
such as ""^^C, ^^O, and ^^Ne are slightly different from the results with the use of 
the reaction network of 464 nucleP^ originating from the differences in the initial 
abundances and adopted reaction rates. The produced mass fractions of elements 
of ^ > 94, which have not been included in the 464 network, amount to about 
10^®. Therefore, differences of 10~^ in mass fractions may be introduced. Although 
there are some differences in the main composition as described in Paper I between 
the original stellar evolution model^ and postprocess hydrostatic nucleosynthesis 
calculations with the reaction network of 464 or 1714 nuclei, the differences are not 
so large for our purpose. 

2.2. S- and p-processes 

Let us focus on the production of elements heavier than iron and the weak s- 
and p-processes in the hydrostatic nucleosynthesis. Figure [T] shows overproduction 
factors X{i)/X{i)Q at the beginning of the core collapse, where X{i) is the value 
averaged over the star including the hydrogen-rich envelope of 38 Mq, and X{i)Q 
is that of the solar system abundances. It is noted that in the estimation of the 
overproduction factors, we adopt the value of 4.55x10^ yr ago for ^*^K, which is a 
long-lived radioactive nucleus (the half-life is 1.25 x 10^ yr) and the present abundance 
is about one order of magnitude less than that when the solar system was born (~ 4.5 
yr ago). There is inconsistency between the initial abundance of ^*^K and the value 
in the estimation of the overproduction factor. However, the overproduction of ^'^K 
is not determined by the initial abundance of ^^K but by the abundance of ^^K as 
described below. Overall, elements of 60 < j4 < 90 are highly overproduced relative 
to the solar ones; the overproduction level ranges over 10^ - 10^, which is similar to 
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Table I. 1714 nuclides contained in the nuclear reaction network for the hydrostatic nucleosynthesis. 



Nuclide 


A 


Nuclide 


A 


Nuclide 


A 


Nuclide 


A 


H 


1 -3 


Cr 


46 -62 


Ag 


100- 


123 


Yb 


162 - 


- 184 


He 


3-6 


Mn 


48-65 


Cd 


103- 


126 


Lu 


165 - 


- 187 


Li 


6-9 


Fe 


50-68 


In 


105 - 


129 


Hf 


168 - 


- 189 


Be 


7-12 


Co 


52-70 


Sn 


108- 


132 


Ta 


171 - 


- 191 


B 


8-14 


Ni 


54-73 


Sb 


111 - 


135 


W 


174- 


- 194 


C 


9-18 


Cu 


57-75 


Te 


113- 


137 


Re 


177- 


- 197 


N 


11 -21 


Zn 


59-78 


I 


117- 


140 


Os 


180- 


-200 





13-23 


Ga 


62-80 


Xe 


120- 


142 


Ir 


183- 


-203 


F 


14-26 


Ge 


64-83 


Cs 


122 - 


145 


Pt 


186- 


-206 


Ne 


17-28 


As 


67-86 


Ba 


125- 


148 


Au 


188^ 


-209 


Na 


19-30 


Se 


69-89 


La 


128- 


150 


Hg 


191- 


-212 


Mg 


21-33 


Br 


72 -91 


Ce 


131 - 


153 


Tl 


194- 


-215 


Al 


23-35 


Kr 


74-93 


Pr 


133- 


156 


Pb 


198- 


-217 


Si 


25-38 


Rb 


76-96 


Nd 


136- 


158 


Bi 


202- 


-220 


P 


27-40 


Sr 


79-98 


Pm 


138- 


160 


Po 


205- 


-222 


S 


29-42 


Y 


81 - 101 


Sm 


141 - 


163 


At 


209- 


-224 


CI 


31-45 


Zr 


83 - 103 


Eu 


143- 


165 


Rn 


212- 


-227 


Ar 


33-48 


Nb 


86 - 106 


Gd 


146 - 


168 


Fr 


215- 


-229 


K 


35-50 


Mo 


89 - 109 


Tb 


148 - 


171 


Ra 


217- 


-232 


Ca 


37-53 


Tc 


91 - 112 


Dy 


151 - 


174 


Ac 


222- 


-234 


So 


39-55 


Ru 


93 - 115 


Ho 


154- 


176 


Th 


225- 


-237 


Ti 


41-57 


Rh 


96 - 117 


Er 


157- 


179 


Pa 


227- 


-239 


V 


43-59 


Pd 


98 - 120 


Tm 


160- 


181 


U 


231 - 


-241 



the weak s-process scenario proposed by many previous studies of the s-process in 
massive stars P'^aJ-ES 

Figure [2] shows the time evolution of mass fractions of selected elements averaged 
over the whole star. ^^Kr, which is one of the representative s-elements, is overpro- 
duced at the helium burning stage. We confirm that the overproduced elements of 
60 < yl < 90 are mainly produced during the helium core burning, and the neutrons 
are mainly supplied by the ^^Ne(a, n)^^Mg reaction (Fig. [2]) as pointed out in pre- 
vious studies (see Ref. [H|) for a recent review). ^^Ne is produced by the sequence 
of ^^N(a,7)i8F(/3-i/)i80(a,7)22Ne reactions. ^^K is produced by the ^'^K{n,-tf^K 
reaction in the helium burning stage. The solar system abundances of ^^K and '^''K 
are 3516 and 0.44 (the present value) (normalized as the abundance of silicon to be 
10^) p3 respectively. Therefore, even if the small amount of ^^K is converted to ^'^K, 
we regard it to be much overproduced relative to the solar value. ^^'^Ta is also over- 
produced in the helium burning stage by the sequence of ^^^Hf(/3~)^^^Ta(n, 7)^'^'^Ta 
reactions (see Fig. [5]). It is noted that the reaction channel of ^^^Hf(/3~)^^^Ta is 
closed if we do not use the /3-decay rates of TY87 because ^"^^Hf is stable in the lab- 
oratory. However, overproduced ^^^Ta decays from a thermally populated ground 
state at the end of the helium burning, and the overproduction level returns almost 
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to the initial value (Fig. [2]). In the helium burning stage, ^^Zr is overproduced more 
than 10 times (Fig.[T]) relative to the solar value by neutron captures, which are faster 
than /3~-decays of ^^Zr (ti/2 ~ 64 d). However, Sr and Y are not so overproduced. 
Although the overproduction levels of the elements of 10 < ^ < 90 do not increase 
significantly after the helium burning, some elements of j4 > 90 are produced after 
that. 

In the carbon and neon burning stages, neutrons are mainly supplied by re- 
actions of i2c(i2c,n)23Mg, ^^C{a,ny^O, ^'^0{-f,ny^O, and ^^Mg{a,ny^Sl The 
^■^Ne(a, n)2^Mg reaction is also activated in the helium burning shell. Since ^^^Hf 
is produced by neutron captures, ^^°Ta is overproduced again relative to the solar 
value via the sequence of ^'^^Hf(/3~)^'^^Ta(n,7)^^'^Ta. 

During the oxygen burning stage, the temperature becomes high (Tc ~ 2x10^ 
K, where Tc is the temperature at the center) and photodisintegrations of s-elements 
activate the p-process in the oxygen- and neon-rich layers. The possibility of the 
p-process in hydrostatic evolution of massive stars was proposed in previous stud- 
jggJ66J,l6i|i,|6BJ.lM|i "Wg find that seed s-elements that have larger mass numbers tend to 
be disintegrated into p-elements by (7, n) reactions, and light p-elements of ^^Se, 
^^Kr, and ^^Sr are produced at the oxygen- and neon-rich layers by (7, n) and 
(7,p) reactions (Fig. [2D. It is noted that ^'^'^Ta is produced more and more by 
the ^^^Ta(7, n)^'^''Ta reaction. Since the solar abundances of ^^''Ta are much smaller 
than those of ^^^Ta, conversion of small amounts of ^'^^Ta leads to the overproduction 
of ^^''Ta relative to the solar value. 

In the later stages, produced p-elements of ^ > 90 are disintegrated by subse- 
quent (7, n) reactions and /3+-decays except for ^^'^Ta. After all, p-elements whose 
overproduction factors are greater than 10 are only Se, Kr, and ^^°Ta at the 
beginning of the collapse (see Fig. [3|). It is noted that (n, 7) and (7, n) reactions are 
in thermal equilibrium at this stage. 

The quantitative assessment of the adopted beta decay rates (TY87) should 
require more experimental data.^^ Moreover, the table of the /3-decay rates of TY87 
covers only the temperature and density of the helium burning stage as mentioned in 
§2.1. Therefore, we have also calculated the nucleosynthesis using laboratory /3-decay 
rates instead of those of TY87. Although overall overproduction levels change only 
by some factors as seen in Fig. HI those of more neutron-rich isotopes of 60 < A < 85 
tend to be decreased. We find that the production of elements for 85 < A < 100 is 
very sensitive to the differences in /3-decay rates, and Sr, Y, Zr, and Mo tend to be 
overproduced. Although the channel of ^'^^]ii{f]~y^'^Ta{n, 'y)^^^Ta is closed without 
TY87, ^^'^Ta is produced by the ^^^Ta(7, n)^^''Ta reaction after the carbon burning 
and remains through (n, 7) ^ (7, n) equilibrium in the later stages. 

Let us summarize the hydrostatic nucleosynthesis. 1) We confirm the weak s- 
process scenario: s-elements of 60 < A < 90 are highly overproduced relative to 
the solar abundances. 2) High overproductions of ^^'^Ta could be attributed to the 
higher density and temperature evolutional tracks of the 70 Mq star. 3) The p- 
process in the oxygen- and neon-rich layers occurs after the carbon burning stage. 
However, the produced p-elements do not remain due to subsequent (7, n) reactions 
and /3"*"-decays. We suggest that for smaller massive stars, the overproductions 
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of ^^"^Ta are significantly decreased because of the lower density and temperature 
evolutional paths. On the other hand, the p-process could become important because 
the produced p-elements of larger mass numbers may survive in low-density and low- 
temperature environments. 



X 



X 

o 




100 120 140 160 180 200 220 240 
Mass number 



Fig. 1. Overproduction factors X{i)/X[i) q against the mass number A averaged over the star 
including the hydrogen envelope at the beginning of the core collapse. Distinguished symbols 
connected by lines indicate the isotopes. 



§3. Explosive nucleosynthesis in a magnetohydrodynamical jet 



In Paper I, we investigated the explosive nucleosynthesis with the nuclear reac- 
tion network including 464 nuclei (up to ^^j^p^JMIi j^ the present paper, we recalculate 
the explosive nucleosynthesis with a much larger reaction network and focus on the 
production of heavier elements, that is, p- and r-elements. In this section, we present 
the review of the explosion model, the method of explosive nucleosynthesis, and the 
results combined with that of the hydrostatic nucleosynthesis. 

3.1. Supernova explosion model 

In Paper I, we have constructed supernova explosion models using a collapsar 
model. ^^'' Here, we briefly summarize the explosion model for the nucleosynthesis 
calculation. We have performed two-dimensional MHD simulations of the collapsar 
model using a nonrelativistic MHD code, ZEUS-2D,l53i'ESl' which is modified^-' for 
handling supernova simulations with a realistic nuclear equation of state (EOS) based 
on the relativistic mean field theory.'^^J 

For a low-density region oi p < 10^ g cm~^, another EOS is connected^ 
smoothly at the density boundary, which consists of the nonrelativistic ions, par- 
tially degenerate relativistic electrons, and radiation. 

We have taken into account neutrino cooling by electron-positron (e ) pair cap- 
tures on nucleons, e^ pair annihilation, and nucleon-nucleon bremsstrahlung. We 
neglect to include both detailed neutrino transport and heating processes, because 
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Fig. 2. Changes in mass fractions of selected elements averaged over the helium core against re- 
maining time before the core collapse (icoii. = 4.22x10^ yr). Upward (downward) arrows denote 
the time of the start (end) of the burning stages. 



the maximuni density remains less than 10 '^ g cm"'^ in our calculations. We discuss 
the effects of neutrino absorptions on the nucleosynthesis in §3.2. BH was mimicked 
as an inner free absorption boundary and gravitational point source with pseudo 
Newtonian potential.'^ We adopted the spherical coordinate (r,0,</>), and the com- 
putation domain was taken from the inner boundaries rjn = 50 - 200 km to 3x10^ 
km, which covers inner oxygen-rich layers. 

The initial presupernova model is the 32 Mq helium core corresponding to an 
-^ms = 70 Mq starp^i which is the same stellar evolution model obtained in the 
previous section. The initial configuration of angular velocity and magnetic field was 
implemented by analytical form with parameters as in the previous study.'^''^''^''^ 
The initial angular velocity is written as follows: 



n (r) = I?o 



r^ -j- r. 



2' 



(3-1) 



where r is the radius from the center, and Hq and tq are model parameters. The 
initial toroidal magnetic field is given in proportional to the angular velocity distri- 
bution as 



S^ {r) = Bo 



^2' 



(3-2) 



where Bq is a model parameter. We adopt the R51 model as in Paper I for the 
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Fig. 3. Same as Fig. [T] but for 35 p-elements. 
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Fig. 4. Same as Fig. [l] but for the results with use of laboratory /?-decay rates instead of those of 
TY87. 



explosive nucleosynthesis, which has the largest amount of the mass end energy ejec- 
tion rates among the investigated models. The parameters of the initial angular 
velocity and magnetic field are Qq = "a s~^, tq = 1500 km, Bq = 5.7 x 10^^ G, and 
Bz = 5 X 10^^ G, where Bz is the initial uniform poloidal magnetic field along the 
rotational axis. The specified parameters of the rotation and magnetic field cor- 
respond to the model of the rapid rotation and the strong magnetic field. Recent 
stellar evolution models indicate^ that if the magnetic field is taken into account, 
the resultant specific angular momentum of the central region becomes smaller than 
that required for the typical collapsar model (j ~ 10^"^ cm'^ s~^ I2D), j^ our simula- 
tion, the jet is triggered by the central magnetic pressure, which grows due to the 
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compression and winding-up effects of the magnetic field, and the amplified mag- 
netic field reaches around ~ 10^^ G. If the magnetorotational instability (MRI}^ 
successfully operates in the core-collapse phase, the magnetic field could be quickly 
amplified to the same level from an initial magnetic field weaker than that ascribed 
in the present paper. However, resolving MRI in a global simulation is very hard 
and not feasible in the present calculations. Therefore, we assume that some mech- 
anisms such as MRI amplify the magnetic field rapidly from a weak initial magnetic 
field and we mimic the situation by simply imposing a strong initial magnetic field. 
Note that the reached magnetic field strength ~ 10^^ G is comparable to that at 
saturation due to MRL^B 

The resulting total ejection mass and explosion energy at the end of the sim- 
ulation {tf = 1.504 s) are 0.124 Mq and 3.02x10^*^ erg, respectively. The specific 
angular momentum j after the formation of a disklike structure is about 5 x 10^® 
cm^ s^^ at the radius of 500 km. The disk extends to about 1000 km from the cen- 
ter at the accretion phase. After the strong jet formation, an expanding bow shock 
is generated at the outer region of the disk and the matter residing in the region 
begins to expand outward along the equatorial axis (see the description of Fig. [5] in 
§3.2). However, even after the formation of the jet and the expanding bow shock, 
the accretion continues at the inner edge of the disk and the accretion rate main- 
tains the value of 0.1 - 1 Mq s~^ with a few factors declined in 1 s. Therefore, our 
model can be regarded as a collapsar model. It is noted that the jet obtained by the 
simulation is mildly relativistic (< 0.1 c) and baryon rich (> 0.1 Mq). In contrast, 
ultrarelativistic GRB jets should be baryon poor (~ 10~^ M(^^). Moreover, the 
event rate of mildly relativistic jets could be larger than those of normal GRB jets.l^ 
Therefore, both the ejected mass and event rate suggest that the contribution to the 
chemical evolution of galaxies could be large compared with that of ultrarelativistic 
ones. 

3.2. Computational method of nucleosynthesis inside the jet 

For calculating the nucleosynthesis inside the MHD jet, 20,000 tracer particles 
are distributed over the computational domain between 1000 km and 3 x 10^ km 
from the center, which covers initially the region from around the iron core surface 
to the inner oxygen-rich layer. All the matter initially located at radii smaller than 
1000 km is absorbed into the inner boundary. The Lagrange evolution of density 
and temperature of each tracer particle can be obtained from the method described 
in Refs.lHO]) and[T]), by which we calculate the nucleosynthesis and follow the change 
in composition. Figure O shows the distribution of the tracer particles at the end 
of the simulation (tf = 1.504 s). The particles initially located at the inner iron 
core. Si-rich layers, and oxygen-rich layers are indicated in red, green, and blue, 
respectively. Some fractions of the particles initially located at the inner iron core 
are ejected by the jet. It should be noted that in Fig. [5l we can see a "blank" 
region in which there are no particles in the equatorial regions. The blank region 
corresponds to the expanding region after the jet formation as mentioned in §3.1. 
The density of the equatorial blank region is p < 10^ g cm^^. In the tracer particle 
method, the distributions of the particles tend to be sparse in expanding and low- 
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density regions. We can also see relatively sparse regions in polar regions. The 
blank is just the problem of the tracer particle method and the region does not 
affect the nucleosynthesis outcome in the ejecta. Additionally, we have performed a 
convergence test of the nucleosynthesis results inside the jet by changing the number 
of distributed particles and confirmed that the results are not changed much by 
the difference in Paper I. Particles that appear deep inside the original iron core 
go through high-density and high-temperature regions; if the temperature is greater 
than 10^° K, the compositions are determined under nuclear statistical equilibrium 
(NSE) condition: they are obtained from the values of (p, T,Y(,). Since these particles 
suffer from electron captures, we need to calculate the change in Y(. of the ejected 
tracer particles due to the weak interactions of e^ captures and /3^ decays until the 
last stage of NSE before the network calculation. The change in Y^. is given bj^^ 

^ = 5;(A+-A„)y„ (3-3) 

i 

where A+ represents the j3~ and positron capture rates and A_ represents the /?"*" 
and electron capture rates. Figure [6] shows the distribution of the ejected mass in 
Mq against the electron fraction of ejected particles at the end of NSE (l"e,/)- After 
the end of NSE, the nucleosynthesis calculation is done along the Lagrange evolu- 
tion of each particle by using a large nuclear reaction network. Since the time of 
hydrodynamical simulation is insufficient to follow the decays of radioactive nuclei, 
the density and temperature profiles of particles are extrapolated assuming adiabatic 
expansion as in Refs.[T|) and l56p . We continue the nucleosynthesis calculation of the 
radioactive decays until ~ 10^*^ yr after the explosion. Note that we neglect the 
feedback of energy generations due to nuclear processes such as photodisintegrations 
in the hydrodynamical calculation because our nucleosynthesis calculations are just 
postprocessing. The effects of neutrino absorptions on the evolution of Yg should be 
noted here because we neglect the effects in the nucleosynthesis calculations, but it 
could be critical for the nucleosynthesis outcome. The inner edge of the accretion 
disk of the R51 model has p ~ 10^ g cm^'^ and the estimated neutrino luminosity 
ranges from ~ 10^^ to ~ 10^^ erg s~^, which is about one order of magnitude smaller 
than that of canonical core-collapse supernovae. Therefore, the neutrino absorptions 
could not be effective. Note that the range of the neutrino luminosity does not 
differ much from that of other models in Paper I. Additionally, although the progen- 
itor and specified initial angular momentum and magnetic field distributions differ 
from that of our model, recent nucleosynthesis calculations in a magnetorotation- 
ally driven core-collapse supernova model with an approximate neutrino transport 
scheme including effects of neutrino absorptions on Yg have revealed^' that the peak 
distribution of Yg in the ejecta is shifted from ~ 0.17 to ~ 0.15 between with and 
without neutrino absorptions, but the results of the nucleosynthesis are not much 
affected by the difference. From the above considerations, the effects of neutrino 
absorptions on Y^ in our model would be small and we neglect the effects. 

To investigate the heavy-element nucleosynthesis including the p- and r-processes, 
we calculate the nucleosynthesis along some Lagrange tracks of the ejected par- 
ticles explained above using the large nuclear reaction networUS^i including 4463 
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at the end of NSE stage. 



nuclei (up to ^^^Am), where the reaction rates are mainly based on the REACLIB 
database^ 'l^ and the adopted theoretical mass formula is the extended Thomas- 
Fermi plus Strutinsky integral (ETFSI).'^ It is noted that the reaction rates in the 
hydrostatic nucleosynthesis calculations are based on the JINA REACLIB database, 
which use the theoretical mass formula of the finite-range droplet model (FRDM) P^ 
Therefore, there could be some inconsistency between the hydrostatic and explosive 
nucleosynthesis calculations. However, taking into account the marked uncertainty 
of the reaction rates far from the valley of the nuclear stability, the differences are 
acceptable in the present purpose. 

3.3. P- and r-processes 

To investigate the total yield of the ejecta, we combine the results of the hy- 
drostatic and explosive nucleosyntheses. We assume that all the unshocked matter 
located at larger radii than the jet front in the star (including the hydrogen envelope) 
oi 9 < 15° is successfully ejected by the jet, keeping the precollapse compositions 
unchanged. 

As shown in Fig. [6l the mass distribution of the ejecta tends to decrease as 
the electron fraction l^j decreases. The lowest value of Yg,/ is 0.192 among the 
ejected particles. We find that the low-Y^j particles are ejected from deep inside the 
disk, which are originally located from the edge of the iron core to the inner Si-rich 
layer and fall into near the BH due to gravitational collapse. The low- Ye,/ particles 
strongly suffer from electron captures, which reduces the electron fractions of the 
particles. The distributions of compositions of the particles are initialized by the 
final state in NSE. On the other hand, the particles that do not suffer from nuclear 
burning are just pushed up by the inner jet at larger radii. Such particles maintain 
the precollapse compositions. Therefore, the precollapse abundances are crucial in 
part to determine the total compositions of the ejected matter. 



14 M. Ono et al. 

The final overproduction factors X{i)/ X[i) q averaged over the ejecta against 
the mass number are shown in Fig. [71 Symbols connected by lines indicate the 
isotopes. Note that neutron-rich elements of 45 < ^ < 55 and 60 < A < 160 are 
highly overproduced relative to the solar values. Figure [8] shows the abundances of 
ejected particles that have different electron fractions at the end of the NSE stage. 
The overproduced elements of 140 < A < 200 originate from the ejected matter, 
which has a lower Y^j of around 0.2 (Fig. [8]), where the ejected materials undergo 
r-process nucleosynthesis. On the other hand, the ejected particles of Y^j ~ 0.3 
produce elements of 60 < j4 < 90. The overproduced elements of ^4 > 90 are 
primarily synthesized in the jet except for ^^''Ta. The overproduction factors have 
a peak aX, A = 195 (the neutron magic number of 126). It should be noted that 
our jet model cannot considerably produce the elements around the third peak. In 
contrast, they are significantly produced in the study by Fujimoto et al.,'^''^ which 
is attributed to the different distribution of the lowest part of Y^j of ejecta. In 
Fujimoto et al., the particles with Y^j ~ 0.1 are also ejected from their collapsar 
model of 40 Mq (see e.g.. Fig. 5 in Ref. [2j)). Since the lowest Yf.j is about 0.2 
in our model, strong r-process and fissions do not proceed. The difference in the 
distribution of Ygj may be ascribed to those of the progenitors and implemented 
initial distributions of the angular momentum and the magnetic field. In particular, 
we include the initial toroidal magnetic field, which may inhibit matter to fall deep 
inside the core and suffer from strong electron captures. 

5-elements of 60 < A < 90 are overproduced significantly, which is due to the 
weak s-process in the hydrostatic evolutional stage. In contrast, no s-elements of 
A > 90 are overproduced, which might be compensated by the products of the main 
s-process in the relatively low mass AGB stars. 

The overproduction factors of 35 p-elements are shown in Fig. [9j P-elements 
whose overproduction factors are greater than 10 are ^^Se, ^^Kr, ^'^Sr, ^^Mo, ^^'^Ta, 
is°W, ^s^Os, i9°Pt, and ^^^Hg. On the other hand, p-elements of ^"^Se, ^^^Kr, ^^^Sr, 
180 W, and ^^''Ta are overproduced in the hydrostatic nucleosynthesis (see Fig. [3|); 
those of ^^Mo, ^^'*0s, ^^"^Pt, and ^^^Hg are mainly overproduced in explosive nucle- 
osynthesis. Underproduced ^^elements in the previous study of the p-process in Type 
II supernovaei^ such as ^^Mo and ^^'^®Ru are not produced in our explosion model. 
Rayet et al.'^''^ have investigated the p-process (gamma process) in parametrized 
supernova explosion models. They concluded that the production of ;^elements is 
very sensitive to the maximum temperature reached during the explosion, and the 
light (iV < 50), intermediate (50 < iV < 82) and heavy {N > 82) p-elements 
are produced in the peak temperature classified as Tg^max ^ 3, 3 > Tg^max ^ 2.7, 
2.5 > Tg^max) respectively, where Tg = T/(10^K). The mass distribution against 
the peak temperatures of the ejecta initially located in the oxygen- and neon-rich 
layers is shown in Fig. [TUl The major part of the ejected matter has Tg^max ^ 2.5. 
Thereby, heavy p-elements are considerably produced in the explosive nucleosynthe- 
sis (Fig. [9|). It is emphasized that light p-elements of ^^Se, '^^Kr, ^^Sr, and ^^Mo are 
primarily synthesized in the ejecta of Tg^max ~ 16 and l^j ~ 0.48 initially located 
in Si-rich layers. Since the peak temperature is very high, protons, neutrons and 
alpha particles are significantly produced by photodisintegrations of heavy nuclei. 
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and the light p-eleinents are produced by a sequence of neutron captures and sub- 
sequent proton captures as described in Fujimoto et al.P the scenario of which was 
originally proposed by Howard et al.!^ Since light p-elements are already produced 
in the hydrostatic nucleosynthesis, the increments due to the primary process in the 
explosive nucleosynthesis are not prominent except for ^^Mo. Note that in Fujimoto 
et al.P intermediate p-elements such as ^^^In, ^^^Sn, and ^^^La are also produced 
by fission only in the ejecta of Yej ~ 0.1. In our model, all the ejected matter has 
Yf,j > 0.2; therefore, fission reactions are not effective and intermediate p-elements 
cannot be produced. 
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Fig. 7. Same as Fig. [T] but for the ejecta with the use of the larger network. 
§4. Summary and discussion 



We have investigated the nucleosynthesis in a massive star of 32 Mq helium 
core with solar metallicity during the stellar evolution and the jetlike supernova 
explosion. In this section, we summarize the results and discuss the uncertainties 
of the production of elements and contribution to the chemical evolution of galaxies 
(§4.1). We also give additional discussion by comparing with observations in §4.2. 

Hydrostatic nucleosynthesis: 1) S'-elements of 60 < ^ < 90 are highly over- 
produced relative to the solar abundances, which is similar to the weak s-process 
scenario proposed in previous studies .'^■'^''^ 2) Although photodisintegrations of 
seed s-elements during oxygen burning produce p-elements, the produced elements 
are disintegrated in the later stages except for ^^"^Ta. 3) Three elements, Sr, Y, and 
Zr, are not much overproduced compared with the solar values except for ^^Zr. 

Explosive nucleosynthesis: 4) Elements of 90 < A < 160 are significantly over- 
produced relative to the solar values. 5) The overproduced elements of 140 < A < 
200 originate from the ejected matter with lower Yg,/ around 0.2, which results in the 
r-process. 6) The p-process produces mainly heavy p-elements {N > 50) because the 
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Fig. 8. Abundances of ejected particles with Ye,/ that have different electron fractions at the end 
of NSE stage. 
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Fig. 9. Same as Fig. [7| but 35 p-elements. 



peak temperatures in the oxygen- and neon-rich layers are relatively low. 7) Light 
^elements are produced as primary ones in the ejected matter, which has a high 
peak temperature. 8) Compared with the previous study of r-process nucleosynthe- 
sis calculations in a collapsar model of 40 Mq by Fujimoto et al.ji''''^ our jet model 
cannot considerably produce both the elements around the third peak of the solar 
r-elements and intermediate p-elements. This may be attributed to the differences 
in the progenitor and the specified initial angular momentum and magnetic field 
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Fig. 10. Ejected mass originated from initial oxygen- and neon-rich layers against the peak tem- 
perature r9,n,ax (Tg = T/ (lO^K)). 



distributions. 

After all, our supernova explosion model results in neutron-rich elements of 70 < 
A < 140 and weak s-elements of 60 < j4 < 90. The origin of other underproduced 
elements would be ascribed to different explosion mechanisms of supernovae. 

Here, we try to deduce the qualitative constraint on the event rate for our 
explosion model. Let us assume canonical supernova explosion to be spherical and/or 
neutrino-driven aspherical ones such as in Refs. [53|) and I55p. The overproduction 
levels for the elements of 60 < j4 < 160 are 1-2 orders of magnitude higher than 
that for those of 20 < ^4 < 60 compared with canonical ones. If we neglect the 
mass loss, the total ejected mass of our jetlike explosion model is around 2 Mq, 
which is about one order of magnitude less than that of canonical ones. Unless the 
event rate of our model is comparable to or one order of magnitude less than that of 
canonical ones, the elements of 60 < ^ < 160 are too produced to explain the solar 
system abundance pattern. Therefore, the event rate of our model could be one order 
of magnitude less than that of canonical ones. The origin of other underproduced 
elements would be ascribed to different types of supernovae. However, the constraint 
speculated here is not strict. 

4.1. Uncertainties of production of elements related to the chemical evolution of 
galaxies 

We discuss the uncertainties concerning overproductions. To investigate the 
weak s-process in the hydrostatic nucleosynthesis, we adopt /3-decay rates of TY87. 
However, the table of the rates for p and T only covers the temperature and density 
ranges in the helium burning stage. As noted in §2.2, we find that productions of 
neutron-rich elements of 85 < j4 < 100, in particular, Sr, Y, Z, and Mo, are very 
sensitive to the adopted /3-decay rates. Therefore, it is urgent to construct a table 
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of the /3-decay rates, which covers all evolutional stages. 

In the explosive nucleosynthesis calculation, we use the theoretical reaction rates 
based on ETFSI mass formula for elements far from the valley of the nuclear stability, 
in which no experimental cross sections are available. The produced abundance 
pattern should depend on an adopted mass formula.'^ 

We assume that all the matter that has larger radii than the jet front in the 
star oi 9 < 15° is successfully ejected by the jet. However, the angle from which the 
matter should be ejected is rather uncertain. Although we neglect the effects of the 
mass loss, a star of 70 Mq with solar metallicity should suffer from a significant mass 
loss.'SSl' If the mass loss is effective, almost all the hydrogen-rich envelope should be 
ejected as the stellar wind. If the above uncertainties are taken into account, the 
overproduction levels could become one order of magnitude less than those described 
in the previous sections. However, the overall abundance patterns do not change 
qualitatively except for the overproduction levels. 

To investigate the possible contribution of our explosion model to the chemical 
evolution of galaxies, we estimate overproduction factors averaged over progenitor 
masses weighted on the basis of Salpeter's stellar initial mass function (IMF).'^ The 
IMF averaged overproduction factor of the element i, /imf (^), is defined as 

Ma 

X {i, M) /cj (M) M (f) (M) dM / X{i)Q 

/IMF (0 ^ '''' M, ' (4-1) 

V / X (j, M) /ej (M) M cj) (M) dM 

where X{i, M) is the mass fraction of the element i in the ejecta, M the initial 
progenitor mass, /cj the mass ratio of the ejecta to the initial mass and (M) oc 
j^^-2.35 ^]^g number of stars within the mass range between M and M + dM, which 
is Salpeter's IMF. Let us adopt spherical postexplosion models from Rauscher et 
al.lSB for progenitors of 15, 19, 20, 21, 25, 30, 35, and 40 M^ with the solar 
metallicity and take our model for 70 Mq. During the integrations in equation 
(j4-ip . necessary values of X(i,M) and /ej are obtained by interpolation of sample 
values. We take Mi to be 15 Mq and M2 to be 70 Mq. For our 70 Mq model, we 
assume all of the hydrogen envelope to be ejected as the stellar wind, because 70 
Mq with the solar metallicity may strongly suffer from the mass loss .'SSI Figure [11] 
shows the IMF averaged overproduction factors. We also show the overproduction 
factors averaged only from 15 Mq to 40 Mq in Fig. [12] for reference. We can see 
relatively high overproduction of ^'^K relative to the solar value, which arises from 
the large enhancement of the star of 20 Mq (see Fig. 4 in Ref. [53])). We also 
recognize relatively larger overproduction factors for elements of 60 < A < 90. The 
overproduction level of 60 < vl < 90 is slightly enhanced by a few factors due to 
the inclusion of our 70 Mq model (Fig. [TT]) compared with that averaged over the 



*' We take the nucleosynthesis data of postsupernova models from the web site: 
EttpT77homepages.spa.uinn.edu/~alex/nucleosynthesis/RHHW02.shtml Note that the data of 
30, 35, and 40 Mq are not yet published, and necessary data for the integration in Eq. H4-l|) are 
obtained by interpolation. 
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range only between 15 AIq and 40 Mq (Fig. [T2]) . Recall that the production of 
elements of 60 < A < 90 results from the weak s-process. Overproduced neutron- 
rich elements of 90 < A < 160 are attributed to the 70 Mq model in which the 
elements are primarily synthesized in the explosive nucleosynthesis by the r-process. 
P-elements of 110 < A < 200 are overproduced to the solar values, which mainly 
come from the p- process in 15-40 Mq stars. ^'^'^Ta is highly overproduced in our 
70 Mq model. Therefore, it is interesting whether our model could contribute to the 
solar ^^''Ta abundance. If we average the overproduction factors over the mass range 
only between 15 Mq and 40 Mq (Fig. [T2]) . the IMF averaged overproduction factor 
of ^^'^Ta is 1.18 in logarithmic scale. On the other hand, if we include the 70 Mq 
model in the integration in Eq. (j4-l|) . the overproduction factor becomes 1.23, which 
corresponds to a 12% increase. Therefore, the contribution of 40 - 70 Mq to the solar 
^^^Ta abundance is negligible in view of uncertainties of the present study. Overall, 
our model contributes to the solar weak s-elements of 70 < A < 90 and neutron- 
rich elements of 90 < ^ < 160. However, we should treat the results with caution 
because our jetlike explosion model is only one specific set of parameters of angular 
momentum and magnetic field distributions, and the fraction of such an aspherical 
explosion is highly uncertain. As suggested by Fujimoto et al.,'^ the r-process does 
not occur in a less energetic jet and the jet properties depend on specified parameters 
of the initial angular momentum and magnetic field distributions. If the fraction of 
jetlike explosions among the progenitors is less than 10%, the contribution of jet- 
induced nucleosynthesis above 40 Mq to the chemical evolution of galaxies would 
be minor. Therefore, we should regard the contributions of our model deduced here 
as the upper limits. In addition, we investigate only the progenitors with solar 
metallicity. As a consequence, some simulations of the chemical evolution of galaxies 
are required to ascertain the contribution especially for elements whose productions 
depend on the metallicity, which is beyond the scope of this paper. 
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Fig. 11. Same as Fig. [T] but for averaged over progenitor masses weighed on the basis of Salpeter's 
IMF with the use of spherical explosion models from Rauscher et al. for 15, 19, 20, 21, 25, 30, 
35, and 40 Mq progenitors and our explosion model for 70 Mq. 
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Fig. 12. Same as Fig.[Tl]but for only from 15 to 40 Mq. 



4.2. Comparison with abundances in metal-poor stars and the possibility for LEPP 

Although our progenitor is assumed to have the solar metallicity, the production 
of r-elements, which are primarily synthesized without seeds, does not depend on the 
metallicity. Therefore, it is worthwhile to compare our results with the abundances 
of extremely metal-poor stars, which are not affected seriously by the s-process. It 
is noted that the overproduced elements of Z > 38 are primarily produced in the 
explosive nucleosynthesis. Abundance ratios relative to Sr against the atomic number 
are shown in Fig. [131 Solid and dashed lines denote abundance ratios of our model 
and that of solar system r-process elements,'^ respectively. Circles, triangles, and 
squares indicate the values of observations of very metal-poor stars CS22892-052,^ 
HDSBGOQSSj' and HD122563,'^ respectively. The symbols without error bars indicate 
the upper limits. We can see that the abundance pattern of CS22892-052, which is 
a typical r-process-rich ([Eu/Fe] > 1) star, coincides well with that of the solar r- 
element pattern, although some exceptions are recognized. The r-process-poor stars 
([Eu/Fe] < 1), HD88609 and HD122563, have a clearly decreasing trend as the atomic 
number increases. It is noted that although [Eu/Fe] values of r-process-poor stars 
are low compared with those of the r-process-rich ones, abundances in HD88609 

([Fe/H] 3.0) and HD122563 ([Fe/H] ~ -2.7) should come from the weak r- 

process because the sources of the s-process, AGB stars, have not had sufficient time 
to evolve before the formation of such metal-poor halo stars .'221' The abundances of 
the ejecta of our model show a decreasing trend in proportion to the decrease in the 
atomic number, which is similar to the case of HD88609 and HD122563.'^ This 
result can be attributed to the decrease in the ejected masses at lower values of Y^j 

(Fig. ED. 

The abundance ratios relative to the solar values are summarized in Table HIl 
where the observational values are taken from Sneden et aii^ for CS22892-052 
and Honda et alM) for HD88609 and HD122563. The values of [Sr/Fe], [Y/Fe] 
and [Zr/Fe] are similarly observed and the ejecta also has the same tendency. While 
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[Sr/Eu] of the r-process-rich star CS22892-052 ([Sr/Eu] ^ —1) is very small compared 
with that of the r-process-poor stars HD88609 and HD122563 ([Sr/Eu] ~ +0.3), the 
value of the ejecta ([Sr/Eu] ~ —0.4) is closer to that of the r-process-poor stars than 
to that of the r-process-rich stars, which reflects the declining trend of abundances 
as the atomic number increases (Fig. [T3]) . 

We find that Sr-Y-Zr isotopes are primarily synthesized in the explosive nucle- 
osynthesis by a similar process of primary synthesis of light p-elements as described 
in §3.3. The ejected matter of Tg^max ~ 16 and Y^j ~ 0.45 produces most isotopes of 
Sr-Y-Zr, which is more neutron rich than that in the case of primary light p-element 
synthesis. In such high peak temperature and density, there exist protons, neutrons, 
and alpha particles, where a lot of neutrons are produced by electron captures. After 
the temperature decreases to 8x10^ K, a sequence of neutron captures and /3~ -decays 
produces slightly neutron-rich Sr-Y-Zr isotopes from lighter elements. After the ex- 
haustion of neutrons, proton captures and gamma processes follow. Recall that ^^Zr 
is overproduced relative to the solar value in the hydrostatic nucleosynthesis (Fig.[l]). 
88gj.^ 89y^ and 91.92,94, 962];- a,re highly overproduced relative to the solar values due 
to the primary process in the explosive nucleosynthesis. It is noted that ^^Zr is more 
produced in the explosive nucleosynthesis. 

Travaglio et al.^^ have suggested that based on a galactic chemical evolution 
(GCE) model, a primary process from massive stars (LEPP) other than the general 
s- and r-processes is needed to explain 8% of the solar abundance for Sr and 18% 
of the solar Y and Zr abundances. In their GCE model, the yields of the s-process 
have been derived from AGB models and they have also added a small contribution 
(~ 10% of the solar ones) from the weak s-component for Sr. It is emphasized that 
the contribution from the r-process has been deduced from the very r-process-rich 
CS22892-052,'i3 that is, contributions to r-elements from r-process-poor stars like 
our explosion model have not been included. In our explosion model, the ejecta 
has a larger [Sr/Eu] than that of r-process-rich stars, and Sr, Y, and Zr are mainly 
produced by the primary process. Therefore, our explosion model could be one 
of the sites of LEPP. However, the calculations are limited to only one model of 
the progenitor with the solar metallicity and specific set of parameters of initial 
distribution of magnetic field and angular momentum. As suggested in Ref. [2]), 
ejected masses of r-elements depend on the jet properties such as the explosion 
energies. Therefore, the effects of explosion models on the chemical evolution of 
galaxies remain uncertain and should be studied in the future. 

Table II. Abundance ratios relative to the solar values for extremely metal-poor stars and our 
model indicated by "Ejecta". 





[Fe/H] 


[Sr/Fe] 


[Y/Fe] 


[Zr/Fe] 


[Eu/Fe] 


[Sr/Eu] 


CS22892-052 


-3.1 


+0.6 


+0.44 


+0.78 


+1.64 


-1.04 


HD88609 


-3.0 


-0.05 


-0.12 


+0.24 


-0.33 


+0.28 


HD122563 


-2.7 


-0.27 


-0.37 


-0.10 


-0.52 


+0.25 


Ejecta 


- 


+1.50 


+1.79 


+1.65 


+1.90 


-0.4 
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Fig. 13. Abundance ratios relative to Sr: calculated ejecta (solid line), solar system r-process ele- 
ments (dashed line) (Simmerer et al.'^) and three very metal-poor stars.'^''^ Symbols with 
error bars indicate the observations of their abundances. Triangles, squares and circles represent 
HD88609, HD122563, and CS22892-052, respectively. Symbols without error bar indicate the 
upper limits. 



Acknowledgements 

We would like to thank A. Heger and his collaborators for offering their data. M. 
Ono thanks N. Nishimura for stimulating discussion. K. Kotake is grateful to K. Sato 
for continuous encouragement. This work has been supported in part by Grants- 
in-Aid for Scientific Research (Nos. 18540279, 19104006, 20740150, 22540297 and 
24540278) from the Ministry of Education, Culture, Sports, Science and Technology 
of Japan. 



Appendix A 

Thermal Population of Ground and Isomeric States of ^^''Ta 



^^'^Ta is one of the rarest isotopes in the solar system and it has a long-lived 
isomeric state. The isomeric state has J'" = 9~ and the half-life r^^ is 1.2x10^^ 
yr, while the ground state has J'" = 1"*" and r^^ — 8.152 h. We treat specially the 
reaction rates concerning ^^"^Ta by a method similar to that described in Ref. [53]). 

Because of the selection rule for the spin and parity, the isomeric state (^^''^Ta) 
cannot directly decay into the grand state i^^^^Ta). However, if the temperature is 
sufficiently high, ^^"^"Ta can decay into the ground state through thermally excited 
states. In stellar interiors during some burning stages and in supernova explosions, 
the two states of ^'^'^Ta are thermally populated. In thermal equilibrium, the popu- 
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lation ratio Piso of the isomer relative to the ground state is given as^ 

^g-0.8738/T9^ (A-1) 



(2Jiso + 1) exp(-^iso/A;r) _ 19 ^.Q.gyas/Tg 



(2Jgs + 1) 

where Jgs and Jjso are the spins of the ground and isomeric states, respectively, and 
Eiso the excitation energy of the isomer, Tg = T/10^ K. If the temperature decreases 
to a critical temperature Tcrit, the two states no longer interact with each other. In 
the explosive burning scenario, the critical temperature is crucial to determine the 
amount of ^^'^™Ta. Belie et al^ have derived the effective decay rate of ^^'^™Ta 
by photoactivation experiments. We simply assume Tcrit as 0.35x10^ K from the 
temperature-dependent decay rate (Fig. 4 in Ref. 19 1|) ). The effective ^^'^Ta rates of 
neutron-induced reaction and /3-decays are given by 

■^eS ^ /gs^gs + /isc^iso) \-^''^) 

where /gs and /iso are the fractions of the ground and isomeric states, respectively. 
For this effective decay rate, we derive Ags from the half-life of the ground state and 
Aiso from Belie et al.l^ For the neutron capture of ^^'^Ta, Ags is taken from Ref. [53|) 
and Aiso is taken from the JINA REACLIB database.l^ If T < Tcrit, we make all 
the ^^°§Ta decayed by hand and set /gs to be 0. 
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